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Abstract 

We study the hydrodynamic evolution of cosmological background neutrinos. 
By using a spherically symmetric Newtonian hydrodynamic code, we calculate 
the time evolution of the density profiles of neutrino matter in cluster and 
galactic scales. We discuss the possible observational consequences of such 
evolution and the resulting density profiles of the degenerate neutrino 'stars' 
in galaxies and clusters. 
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I. INTRODUCTION 



The observations of neutrino oscillations indicate that at least one flavour of neutrinos 
have non-zero rest mass, and at least one type of sterile neutrinos with mass of eV order 
must be taken into account for the three scales of mass squared differences if the LSND 
experiment is confirmed [1]. A huge number of primordial neutrinos were produced in the 
early universe [2], [3], and they were cooled down as the universe expanded, decoupling from 
baryons much earlier than the photon decoupling epoch. The temperature of the neutrinos 
dropped to below 0.1 meV after several billions of years, and so if their rest mass rrii, is an 
eV or so, the neutrinos have become non-relativistic. Sterile neutrinos may be produced 
in the early universe by the oscillations of the active neutrinos [4]. These sterile neutrinos 
may have rest mass greater than 10 eV and may play a role in structure formation. Also, a 
neutrino rest mass rrii, of the eV to keV order can account for most of the dark matter, which 
dominates the mass in galaxies and clusters [5]. In this article, we study the gravitational 
collapse of such non-relativisitc cosmic background active and sterile neutrinos, and we 
discuss the formation of degenerate neutrino dark halos in cluster and galactic scales. It has 
been proposed that such degenerate neutrino 'stars' can be found in centers of galaxies [6] , 
[7]. 

II. COSMIC BACKGROUND NEUTRINOS 

In standard cosmology, the neutrinos decoupled at around 1 MeV, and the present value 
of their number density rii, is about 100 cm~^ for each type of neutrinos [8]. Therefore, we 
can write the number density of the primordial neutrinos as 

100 -3 /.N 

rii, cm (1) 

where a is the cosmic scale factor, set to be 1 at present. If the neutrinos were ultra- 
relativistic, the neutrino temperature would be 



where Tq 1.95 K is a constant [8]. When the neutrinos cooled down, they became non- 
relativistic and T^, ~ a~^. Suppose at a certain a — ttc and T^, — Tc, kT^ — rrii^c^, and the 
neutrinos become non-relativistic. We have 

rp Tl 

where Ti is the present neutrino temperature. Therefore, we can express in terms of 
m,^(in eV): 

. ^-^K. (4) 

for the non-relativistic regime. We will use this relation to check the degeneracy of the 
neutrinos at the structure formation epoch in the next section. 

In order to estimate the time for the beginning of structure formation, we should consider 
the Jeans mass of the primordial neutrinos, which is defined as 

Att 

Mj = -p^X'j, (5) 

where pi, is the mass density of the neutrinos and Aj is called the Jeans length which is 
defined by 

with e the mean energy density and Pi, the pressure of the neutrinos. When the mass of an 
object is larger than its Jeans mass, the object will collapse gravitationally to form structure. 
The mass density of the neutrinos in the zero Ti, limit is given by 

Pz. = ^ xpi'^xl + l)^l + xj,-ln(xF + ^l + xl^ , (7) 

where X is a constant and xp is defined as 



2Trh 
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Xi? <^ 1 and Xp ^ ^ represent the non-relativistic and ultra-relativistic regimes respectively 
for the neutrinos. By combining Eqs. (5)- (8), we obtain a relation between Mj and a for a 



given rrii, (see Fig. 1). If the Jeans mass is of galactic scale (about lO^^M©), the value of a is 
about 0.0585 for m^, — 10 eV but greater than one for m^, = 1 eV. For cluster scale (about 
IO^^Mq), the corresponding values of a are 6 x 10"*^ and 0.0125 for m,y being 10 eV and 1 
eV respectively. 

By looking at the Jeans mass only, we find that clusters will start to form first and then 
the galaxies, which is known as the top-down scenario. This disagrees with the observation, 
which indicates that the formation of smaller structures occurred earlier in time [5]. How- 
ever, we should compare not only the time that the structures start to form, but also the 
duration of the formation of structures in order to determine whether neutrinos can still 
form structures and agree with the observation. 

We believe that different masses of neutrinos may correspond to different structures- 
the larger the neutrino mass is the smaller the structure and vice versa. The radius of a 
hydrostatic neutrino star is given by [9] 



where M is the total mass of the neutrino star. The corresponding R are Mpc, kpc and 0.01 
pc for nil, — 1 eV, 10 eV and 1 keV respectively if M = lO^^M©. Therefore, we assume that 
neutrinos with rest mass of 10 eV or above will dominate in galaxy scale and those with rest 
mass 1 eV or below will form structures in cluster scale. 

III. HYDRODYNAMICS OF NEUTRINO STAR FORMATION 

The formation of degenerate heavy neutrino stars was discussed in Ref. 10 using the 
analogy with an interacting self-gravitating Bose condensate. Here, we use a hydrodynam- 
ical code to simulate the formation of neutrino stars in galaxies and clusters. We use the 
Lagrangian mass coordinate to write the spherically symmetric hydrodynamic equations 





[11]: 



du _ 2 \dP dQ 



dt dm dm 



Gm{r) 



(10) 



4 
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Here, u is the fluid velocity, p the mass density, m the enclosed mass, and Q is the artiflcial 
viscous stress. The above equations indicate that the neutrinos are under gravitational 
attraction and influence of the pressure P. There exists neutrino degeneracy pressure as 
the neutrinos become degenerate after they have gravitationally collapsed. The criterion of 
degeneracy is given by 

(m kT 

We can use Eq. (4) to approximate T^, and set a lower limit for the mass density in degenerate 
state. The neutrino degeneracy pressure is given by 



471^^2 / 3 \ 2/3 
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In fact, the neutrinos are coUisionless except for the degeneracy pressure. Therefore, we set 
the effective temperature of the neutrinos to be zero because the neutrino temperature is 
much less than their Fermi energy. Here, we consider non-relativistic degenerate neutrinos, 
the Fermi momentum for which is much smaller than their rest mass. In the hydrodynamic 
code, we use an artiflcial viscous stress to smooth out possible shock waves. The artiflcial 
viscous stress is defined as: 



Q = -a 



du 1 (9 , 2 X 
dr 3r^ dr 



(14) 



where a is the scalar artificial viscosity defined by 



(15) 

with / being a constant length which will be defined later. In the hydrodynamic simulation, 
we use a finite difference scheme to evaluate the above equations [12]. We represent all 
the variables as: u — > u^'^^^'^ — > r^^^^"^ ■, p — > Pfe+i/2)-P ~^ -ffc+i/2' where the upper index 
represents time step and the lower index labels mass shell, and we let 



I' = CQ(Ar,+v2)^ (16) 

where Arfe+1/2 is the thickness of the k*"^ mass shell and Cq is an adjustable parameter set to 
be between 2 to 3, to spread the shock front into several zones. The difference between two 
consecutive Lagrangian coordinates is the mass of a zone, Am Amk+1/2 = i^k+i —trik- In 
addition we use the Courant stability condition and test for convergence to make sure that 
the computer code works effectively. 

In our simulations, we set the initial density profiles to be uniform. We use IO^^Mq and 
IO^^Mq as the Jeans mass to find the values of a at which gravitational collapse begins for 
different rrii,. The initial density po can be obtained from the standard cosmology: 

Po = (17) 

where Vt^ and pc are the cosmological density parameter of the neutrinos and critical density 
of the universe respectively. All the parameters used in the simulations are shown in Table 
1. Also, for simplicity, the initial hydrodynamic velocity of the neutrinos is set to be zero. 

In fact, after the epoch of decoupling, the neutrinos were not fully degenerate and had no 
interaction with other neutrinos, baryons, or photons. Therefore, the pressure is very small 
and not enough to prevent gravitational collapse. However, when the neutrinos continue 
to collapse, their density will increase until they become fully degenerate. Further increase 
of the number density will boost up their degeneracy pressure until the pressure gradient 
balances the gravitational attraction and the collapse is stopped. If the neutrino mass is 
too large, the degeneracy pressure will be small and the neutrinos will collapse to a high 
density within a very small region with a high speed. Therefore, the neutrinos will rebound 
back with a high velocity. If the neutrino mass is too small, the neutrinos would remain 
relativistic and the gravitational attraction is not enough to bind them together. In other 
words, there only exists a range of neutrino mass to form a stable neutrino star. In Fig. 2-5, 
we demonstrate three cases with the parameters in Table 1. In Fig. 2, the resulting density 
becomes smaller as time evolves, which means that the neutrinos cannot form structures by 
themselves if is 1 eV. In Fig. 3, we still use = 1 eV but we impose a constant density 
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core (10~^^ gcm~^ within 500 kpc) inside the cluster. The neutrinos can form structure 
with the help of additional mass. In Fig. 4 and 5, we use = 10 eV and 1 keV, and the 
neutrinos can form structures by themselves. The total mass of the structures we used are 
IO^^Mq for m,y — 1 eV and lO^^M© for rrii^—lO eV or 1 keV. These few figures show the 
possibility of neutrinos forming structures in both galaxies and clusters. Also, the time scale 
of formation is smaller as rrii, is larger. The time scales are about 10^^ s, 10^^ s and 10^^ s 
if nil, are 1 keV, 10 eV and 1 eV respectively. Therefore, when one includes the formation 
time, the sequence of neutrino structure formation agrees with the observation, i.e. smaller 
scales form first. 

In Fig. 6, we compare the density profile at hydrostatic equilibrium with the late time 
profile of the dynamical evolution. They appear almost identical except at the outermost 
region, where the difference arises because of some undamped oscillations in the outer part 
of the neutrino star. The density profile is nearly constant at small radius and drops greatly 
at large radius. 

In cluster formation, most of the baryons become hot gas particles and stars in galaxies. 
Their mass is about one-tenth of the total mass of the cluster. We have ignored the effect 
of baryons here in the simulation. Nevertheless, the constant density profile at small radius 
agrees with some observational data of cluster dark matter [13], [14]. 

When the neutrinos begin to collapse by gravitational attraction, they gain kinetic energy 
and the density of the inner part of the neutrino star increases. When the central density 
is high enough, the pressure gradient prevents the neutrinos to collapse further and some 
of the neutrinos begin to rebound back. We find that there are oscillations in the neutrino 
star. In Fig. 3, we can see oscillations of the central density. Also in Fig. 7, we plot the 
density against time at r = 50 kpc with the same parameters as those used in Fig. 4. The 
period of oscillation r is almost constant and the amplitude varies with time. The value r 
in the central region can be found by Taylor expansion of the hydrodynamic equation with 
small displacement, which is roughly given by the following relation [15]: 
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, ■ (18) 

Prom Fig. 7, the period of oscillation is about 1.6 billion years. The density at r = 50 
kpc is about 6 x 10~^^ gcm~^. Therefore, by Eq. (18), r ~ 1.5 billion years, which is 
roughly consistent with the result. Because of the weak coupling among the neutrinos, the 
oscillations probably would not be damped out for a long time. 



IV. DISCUSSION 

We have demonstrated the possibility that the neutrinos can form structures in galaxies 
and clusters. Including the formation time, our model is consistent with the observation but 
not the top-down scenario. Our hydrodynamical simulation results (see Fig. 4) agree with 
recent observed data indicating that the central density of dark matter profile is constant 
[13]. 

Light neutrinos (1 eV) cannot form structure by their own gravitational attraction. How- 
ever they can still form structures under the gravitational field of a heavier component (see 
Fig. 2 and 3). It is well known that light primordial neutrinos (eV order) cannot account 
for all the dark matter as 0.,^ is much less than 0.3 for m,j = 1 eV [16]. There must exist 
heavier particles such as sterile neutrinos which dominate the mass of the dark matter and 
provide such a gravitational field in the clusters. These sterile neutrinos (10 eV order) may 
also form smaller structures in galaxies (about 10-20 kpc), and all these smaller structures 
link up together to form a large structure in clusters. Therefore, neutrinos with rest mass 
10 eV or above are important in structure formation, and we should not neglect their con- 
tribution. The range of possible masses of active neutrinos and sterile neutrinos can both 
be constrained by more data from the distribution of hot gas in clusters and the rotation 
curves of galaxies. 

From Fig. 3 and 7, we can see that the central density of a neutrino star oscillates, 
which may affect the hot gas in clusters. It may force oscillate the hot gas, which may be 



8 



observable. Also, the oscillations in the hot gas may convert the kinetic energy into internal 
energy of the hot gas and thus heat them up, which may at least partially account for the 
heat source of the hot gas. 
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TABLES 



TABLE I. The parameters used in the simulations. 





Po(gcm ^) 


Cut off radius (kpc) 


1 eV 


1.02 X 10-2^ 


539 


10 eV 


9.97 X 10-27 


11.7 


1 keV 


1.41 X 10-^s 


0.022 
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FIG. 1. Jeans mass of primordial neutrinos with mass of 1 eV (solid line) and 10 eV (dashed 
line) respectively vs. cosmic scale factor. Also, indicated are typical cluster and galactic mass 
scales. 
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FIG. 2. Time evolution of the density profile of 1 eV neutrinos, from initially uniform and 




stationary distribution. Total mass = 10 Mq and the time steps are in 10 years. 
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FIG. 3. Time evolution of the density profile of 1 eV neutrinos with a constant density core 
(~ 10~^^ gcm~^ within 500 kpc) of dark matter and total mass 10^^ M0 in the background. The 
time steps are in 10^ years. 
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FIG. 4. Time evolution of the density profile of 10 eV neutrinos. Total mass = IO^^Mq. The 



time steps are in 10^ years. 
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FIG. 5. Time evolution of the density profile of 1 keV neutrinos. Total mass = WP'Mq. The 



time steps are in 10^ years 
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FIG. 6. Comparison between the final density profile of the simulation and hydrostatic density 
profile with the same parameters as those used in Fig. 4. The solid line represents the final density 
profile of the hydrodynamic simulation. The circles represent the density profile of the hydrostatic 
neutrino star. 
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FIG. 7. Oscillation of the mass density of a neutrino star at r = 50 kpc with the same param- 



eters as those used in Fig. 4. 
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